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Self-Consistent Field Theory is applied to a film of cylindrical-forming block copolymer subject 



^ , to a surface field which tends to align the cylinders parallel to electrical plates, and to an external 

j^ , electric field tending to align them perpendicular to the plates. The Maxwell equations and self- 
consistent field equations are solved exactly, numerically, in real space. By comparing the free 

P^ ■ energies of different configurations, we show that for weak surface fields, the phase of cylinders 

o, 

CJ , parallel to the plates makes a direct transition to a phase in which the cylinders are aligned with 

^vq ' the field throughout the sample. For stronger surface fields, there is an intermediate phase in which 

Q^ ■ cylinders in the interior of the film, aligned with the field, terminate near the plates. For surface 

a^: 

CN . fields which favor the minority block, there is a boundary layer of hexagonal symmetry at the 

^^ I plates in which the monomers favored by the surface field occupy a larger area than they would if 

o 



the cylinders extended to the surface. 



I. INTRODUCTION 

With their abihty to self-assemble, block copolymers are a natural choice of material to 
be utilized in the fabrication of devices incorporating periodic arraysi'^'^. A major difficulty 
with the process, however, is that the system rarely forms a single domain of such an array, 
so a desirable long-range order is absent. Furthermore, preferential interactions between the 
blocks and substrate can cause the array to be aligned in a direction other than the one 
desired for use. Alignment of domains can be obtained in several ways, of which the one of 
interest to us is the use of an electric field^. This takes advantage of the fact that the two 
blocks in general have different dielectric constants, so that it is favorable for the array to 
align itself such that there is as little induced polarization charge as possible^. For lamellar- 
and cylinder-forming diblocks, this means that a sufficiently strong field will align the array 
with lamellae or cylinders oriented with their long symmetry axis parallel to the applied 
field, perpendicular to the electric plates. 

If one assumes that the system, in zero external electric field, is characterized by layers 
of cylinders lying parallel to the substrate and, in large external fields, is characterized 
by all cylinders oriented perpendicular to the field, then one would like to know both the 
equilibrium morphology at intermediate field strengths and the phase diagram of the system. 
Two possibilities were suggested by Thurn-Albrecht et ah^. In the first, the electric field 
has essentially little effect on the arrangement of parallel cylinders up to the field strength 
at which the perpendicular morphology becomes the globally stable phase. If this situation 
obtains, one expects that the difference in electrostatic energies of the two configurations 
would be E^Ado, with A the area of each plate, and do the distance between plates. Hence 

— 1/2 

a first-order transition would take place at a critical electric field, Ec, such that Ec (x d^ 
In the second scenario, an intermediate phase is formed, one in which the surface field is 
sufficiently strong that a boundary layer near the plates has a different morphology than 
in the rest of the film. As Thurn-Albrecht et al. reported observing cylinders both in 
parallel and perpendicular orientations over a range of field strengths, the implication was 
that while the cylinders were oriented perpendicular to the substrate in most of the film, the 
boundary layer might consist of cylinders oriented parallel to the surface. The difference in 
electrostatic energies between the perpendicular phase and this intermediate phase is, then, 
not E'^Ado, but E'^At, where t is a length characteristic of the thickness of the boundary 



layer. Presumably it depends only weakly on the film thickness, do- Thus one expects that 
Ec also varies only weakly with thickness for thick films. Thurn-Albrecht et al.— did indeed 
find that the critical electric field, beyond which only the perpendicular phase was observed, 
was relatively constant for large film thicknesses. 

The existence of an intermediate phase is a consequence of the competition between the 
surface fields, which prefer a phase in which the cylinders are parallel to the surface, and 
the electric field which prefers that the cylinders be aligned with it, perpendicular to the 
substrate. If the system were semi-infinite, the intermediate phase would simply be labeled 
a surface phase, one of different symmetry than that of the bulk, a very common situation^. 
Such a phase would be expected to persist in the case of a finite, thick film. However when 
the film is sufficiently thin, the clear distinction between bulk and surface properties can 
no longer be made, and the two orientations compete throughout the film, resulting in the 
elimination of one of them. 

Most theoretical work on block-copolymers in electric fields has focused on lamellar- 
forming phases, as opposed to cylindrical ones, and on the dynamic mechanisms by which 
morphologies could realign, rather than on the equilibrium morphologies themselves. Such 
dynamic mechanisms were the subject of the seminal papers of Amundson et al.-^, and 
of more recent dynamic density functional calculations^iiSiii. The form of the equilibrium 
phases was taken up by Pereira and Williamsi^ and by Tsori and Andelmani^. Both groups 
considered one intermediate phase in which a single lamella at the plates remains parallel 
to them, although surface fields could conceivably lead to a sequence of intermediate phases 
distinguished by the number of lamellae parallel to the surface. Such a discrete increase, or 
decrease, of the number of lamellae with applied field would be analogous to the layer-by- 
layer growth modes observed in adsorbed films^^. Both groups presented phase diagrams 
as functions of applied electric field and film thickness. Although similar, they differ in one 
respect: that of Tsori and Andelman shows the perpendicular phase to be favored always 
at sufficiently large external fields, which is certainly correct. The effects of surface fields 
have also been exploredi2*i^^. As expected, a strong surface field favors the formation of 
an intermediate phase. The evolution of morphology with external field in a single lamella 
has recently been studied by Matsen using self-consistent field theoryA5. 

As for cylindrical-forming diblocks in an external electric field, the only published work 
has again focused on the mechanisms of morphology realignment. This has been studied 



experiment all}*i^, and also theoretically by dynamical self-consistent field tlieoryii*i^. The 
advantage of this technique is that one follows the dynamics of the system, and one can 
clearly see how different are the pathways towards an equilibrium morphology depending 
upon the initial configuration. The actual equilibrium morphologies, which are the focus of 
our interest, are difficult to achieve however, just as they often are in experiment. 

Some insight into the expected phase diagram of the cylinder-forming diblock copolymer 
system can be garnered from the behavior of the lamellar-forming ones. In both systems, 
surface fields prefer an orientation which is different from that preferred by the electric 
field. It should be noted that this holds irrespective of whether the surface fields prefer one 
monomer species. A, or the other, B. Thus we should expect to obtain a phase of cylinders 
parallel to the surface when the electric field is weak, a phase of perpendicular cylinders 
when the electric field is strong, and an intermediate phase for thick films, and ones subject 
to large surface fields. The principal difference between the cylinder- and lamellar-forming 
systems will be in the symmetry of the surface phase. Similarly the principle difference 
between cylinder-formers adsorbed on a substrate that prefers the minority monomer. A, 
and those adsorbed on a surface preferring the majority monomer, B, will be in the specific 
morphology of the surface phase. Otherwise the general form of the phase diagrams of these 
systems is expected to be similar due to the identity of the underlying physics. 

In this paper, we study a cylinder-forming block-copolymer system in external electric 
fields and examine the morphology and phase behavior in the presence of surface fields which 
prefer the minority monomer. In section II, we first show that some general features of the 
surface phase diagram in the electric field-chemical potential plane are easily discerned from 
basic thermodynamics, and without detailed calculation. However thermodynamics does not 
tell us the nature of the intermediate phases, while specific calculation does. We then turn to 
self-consistent field theory, and solve exactly the self-consistent equations and the Maxwell 
equations in the manner used by us previously for a bulk phase^^-. Subsequently, the same 
idea has been employed with a different technical implementation in other system3i^»i^. As 
a consequence of solving the Maxwell equations exactly, we do not assume, as in refs .'^* - ''^'^'^^ , 
that the effect of the electric field is weak; that is, that the fractional variation of the dielectric 
constant throughout the morphology is smalU^. Our results are presented in section III. For 
weak surface fields, or thin films, there is no intermediate phase, and we find the first of the 
scenarios suggested by Thurn-Albrecht et al^. The phase of parallel cylinders transforms 



directly to one of perpendicular cylinders. These cylinders are, of course, affected in their 
structure by the presence of the substrate surface fields, but a cross-section parallel to the 
substrate confirms the hexagonal order, even very close to the substrate. 

For stronger surface fields or thicker films, we find an intermediate phase in which a layer 
adjacent to the plates has a different morphology than that of the field- aligned cylinders 
in the rest of the film. The morphology is not one of cylinders parallel to the plates, but 
has the same hexagonal symmetry as a cross-section of the field-aligned aligned cylinders. 
However the distribution of monomers is altered so that the monomer favored by the surface 
field occupies a larger surface area than it would had the cylinders extended to the surface. 

We conclude with a brief discussion and comparison of our results with experiment. In 
particular we find, in agreement with experiment in the presence of strong surface fields, that 
the phase boundary between the perpendicular phase and the intermediate phase is a weak 
function of film thickness for large thicknesses. Along this locus of first order transitions, the 
perpendicular and intermediate phases coexist. However as we just noted, the boundary layer 
of the intermediate phase does not consist of cylinders parallel to the substrate. Thus parallel 
cylinders do not exist in either phase, and one would not observe in equilibrium, along this 
particular phase boundary, the coexistence of cylinders in both parallel and perpendicular 
orientations. 

II. THEORY 

A. Thermodynamics 

Because we are dealing with a surface film, we consider the excess free energy per unit 
area, or surface free energy per unit area, a, defined as follows. Let Ftot be the total free 
energy of the system of volume Q = d^A such that the free energy per unit volume, 

is the Legendre transform of the energy per unit volume, e^, with respect to the entropy 
density, Sb, number density, ph, and displacement field i.e. 

fb = eb- SbT - pbp - Dfc ■ E, (2) 



with Dft and E the volume-averaged displacement and electric fields. Then the excess free 
energy per unit area, a, is 

^ ^ lim ^kz_I^ (3) 

with the differential 

da = -SsdT - psdfi - DgdEext- (4) 

Here Sg is the excess entropy per unit area, Ps the excess number of particles per unit area, 
and Eext is the externally applied field which is equal to the spatially averaged electric fielcii^. 
The quantity D^ is the magnitude of the surface contribution to the displacement field, or 
surface excess displacement field 

D, = lim ^ /[D(r) - Bt,{r)]d\ (5) 

where Db(r) is the displacement field in the bulk cylindrical phase. The surface excess 
displacement field is in the x direction, normal to the plates. 

From this differential, one easily finds a Claussius-Clapeyron equation pertaining to the 
surface phases. If one plots the phase diagram at fixed temperature in the electric field- 
chemical potential plane, then the slope of the boundary between surface phases, dEext/dfi, 
is given by 

dEext _ ^Ps /^N 

~iir~ ad's' ^^ 

where Aps and ADg are the differences in the excess surface densities and displacement 
fields in the two phases. At zero external field, there are an infinite number of parallel 
phases distinguished by the integer number of layers of cylinders, and a perpendicular phase 
in which the cylinders are perpendicular to the substrate. The transitions between all phases 
are generally first-order; those between the perpendicular and parallel phases are first-order 
due to the difference in symmetry between them; those between different parallel phases are 
generally first order because an additional layer of parallel cylinders cannot be added in a 
continuous manner. As one turns on the electric field, the phase space of the perpendicular 
phase will increase, and that of the parallel phases will decrease. It is clear that the spatially 
averaged polarization will be smaller in the parallel phases than in the perpendicular phase, 
and so the displacement field will decrease discontinuously on entering a parallel phase from 
a perpendicular phase, and will increase discontinuously when leaving a parallel phase. Fur- 
ther, the surface density is presumably a monotonically increasing function of the chemical 



potential. Hence we expect from eq. IHl to find the infinite number of parallel phases to be 
enclosed by phase boundaries which, if no other phases are encountered, will be in the shape 
of wickets. The legs of each wicket meet the E^^t = axis perpendicularly, because ADg = 
in the absence of a field. The top of each wicket is locally horizontal, and the field at which 
it occurs defines the critical voltage at which that parallel phase disappears. In the region 
of phase space in which an intermediate phase appears, the Claussius-Clapeyron equation 
can be applied to the phase boundaries between it and the parallel phases, and between 
it and the perpendicular phase to obtain useful information relating the surface densities 
and displacement fields in the coexisting phases. It can also be applied, of course, to the 
lamellar-forming system a-^^i-^^ . 
Another useful result, 

follows from the differential eq. |3] We note that, whereas that component of the bulk, 
spatially averaged, displacement which is along the direction of the external field must be 
positive^Si, that same component of the excess surface displacement need not be. In fact, 
as we expect that the presence of the surface can only disrupt the perfect alignment of the 
cylinders with the electric field, at least near the surfaces, we anticipate that Ds will be 
negative. It follows from eq. [7|that the excess free energy per unit area will increase with 
increasing electric field. This contrasts with the behavior of the total free energy, Ftot, and 
free energy per unit volume, fb, which decrease with increasing external field. 

B. Self- Consistent Field Theory 

The method we employ has been described previously^ for a bulk system, and simply 
needs known modifications for the case of a surface filnt^i. Thus its presentation here can 
be brief. We consider a melt of n A-B diblock copolymer chains, each of polymerization 
index A^ = A^^ + Nb- The mole fraction of the A- monomers, /^ = N^/N . We also assume 
that the Kuhn lengths of the A and B components are identical, a length denoted a. The 
system is confined between identical plates, a distance do apart, each of area ^, which exert 
a surface field, /i(r), preferentially on one block with respect to the other. The plates are 
normal to the x axis, and the applied field is along the positive direction. The fraction of 
the volume occupied by polymer is denoted $o(r); and is unity in the film except near the 



plates where it falls to zero. 

Self-consistent field theory leads to a free energy, F, which is a functional of unknown 
fields, Wa, Wb, S, and the unknown electric potential, V{r), and a function of temperature, 
T, volume, Q, and area, A. In an ensemble in which the external electric potential, Vq, is 
held HxedM, the free energy can be written 

-— = -\nQ[WA,WB\ 



^ a 



- J dr{xN^A^B - Wa^a - Wb^b 

(VK(r))n , (8) 



fleoK^r) f^^^i ^.2i 



InkBT 

where eo is the vacuum permitivity, A;^ is Boltzmann's constant, and ^^(r) and ^^(r) 
are the local volume fractions of A and B monomers. The function Q[Wa,Wb] is the 
partition function of a single polymer chain subject to the fields W^A(r) and W^s(r), as is given 
below. The field H(r) is a Lagrange multiplier that enforces locally the incompressibility 
constraint, $^(1") + ^siT") = "^o- The dependence on temperature, T, comes from the usual 
Flory interaction parameter, x, which to a good approximation is inversely proportional to 
T. Finally the local dielectric constant is ^(r). A constitutive relation between this local 
dielectric constant and the local volume fractions, $A(r) and $s(r), must be specified. We 
have chosen, as previously^^, that the local dielectric constant be given by its local average. 
Here that choice is expressed 

K{r) = ^olKA^Air) + KB^Bir)] + 1 - $0, (9) 

where ka and kb are the dielectric constants of the pure A and B homopolymer phases 
respectively. Note that near the plates, the dielectric constant approaches unity. We stress 
that, while this choice is clearly correct in the limiting cases of the pure systems and in the 
weak segregation limit, it is simply a choice. While there are theories of the constitutive 
relation appropriate for dilute gases and for critical fluids, little is known of the detailed 
relation for dense, solid, multi-component systems. Hence that of eq. IHl which has the 
virtue of being both reasonable and simple, has been commonly employe d-'^' i ^'^'^i'^^'^i'^^i'^'i'i'^^ . 

The requirement that the free energy functional be an extremum with respect to variation 
of the electric potential, V{r), ^^^(r), WB^r), H(r), and of the volume fractions $A(r) and 
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$B(r) at constant temperature, or x^, leads to the following set of equations: 



= V ■ [eo«:(r)V\/(r)] , 



WB[r 
0A(r 



Qen 



XiV0A(r)+e(r) + /i(r)iV- 

n 6Q 



2nkBT 
2nkBT 



%{v)KA{VV{v)f 

<l>o(r)«:B(VV(r))^ 



Q5wA{r) 

n 6Q 



(10) 

(11) 

(12) 
(13) 
(14) 

(15) 



QSwBir) 

The functions Wa, Wb, S, $^ and $b, which satisfy these equations are denoted by lower 
case letters, wa, wb, ^, (pA, 0b, respectively. The first of these equations is simply Gauss's 
law for a system with no free charge. The free energy within the self-consistent field approx- 
imation, Fgcf, is obtained by substitution of these functions into the free energy of eq. (jHI), 



Fsef (K), T, n, A) = J'iwA, Wb, e, V; T, n, A) 



(16) 



or 
nkBT 



\nQ[wA,WB]-^J dr{xNM^)M^)H{r)M^) + ^^[l-MrmViv)f} . 



(17) 

The nature of the particular system being described by self-consistent field theory is 
specified by the partition function of a representative member of the entities which comprise 
the system, in this case, block copolymers. Therefore QIwajWb] = j d^^ (li^A)/^., where 
g(r, s) satisfies the modified diffusion equation 



and 



^ = ]:Na'V^q-WA{v)q, if 0<s<^ 

OS 6 



^ Na^V^q-WB{Y)q, if /a < s < 1 



(18) 



(19) 



ds 6 
with the initial condition g(r, 0) = 1, and A; is a volume of no consequence here. 

The total density profile, $o(r), and the surface field, /i(r), must still be specified. We 

follow ref.— and choose 



$o(r) 



\[1 — cos(7rx/e)]. 



<a; <e 



and 



= 1, e<x<(io — e 

= |[1 — cos(7r((io — x)/e)], do — e < x < do (20) 

h{r)= ^^^^[l + cos{-Kx/e)], 0<x<e 

= 0, € < X < do — e, 

= ^^^^[l + cos{7T{do-x)/e)], do-e<x<do, (21) 

where e is a small distance given below, and A is the strength of the surface field. A positive 
value of A causes the surface of the system to prefer the A component. 

Because of the three-dimensional nature of the problem, we have chosen to solve the 
self-consistent equations in real-space, following the pseudo-spectral method of Rasmussen 
and co-worker322i2i2i. We have taken a grid which is 192 x 32 x 64 in the x, y, z directions 
respectively. The length e is chosen to be sufficiently small to ensure that it does not affect 
phase behavior. It is discussed further, below. In examining the phase of cylinders parallel 
to the plates, whose axes are in the z direction, we have investigated several arrangements 
of the cylinders, and varied the distances between them, both in the x and y directions, to 
ensure that we have found the minimum free energy configuration. 

III. RESULTS 

The parameter space of our system is large, as its state is specified by five variables, 
the choice of which is not unique. They are the temperature, or equivalently, x^^ the 
applied electric field, Eext, the chemical potential, /i, or the film thickness do, the surface 
field strength, AA^, and the copolymer architecture, Ja- In order to demonstrate the var- 
ious morphologies, we have chosen to fix x^ = I85 ^ind the architecture at Ja = 0.29 
which corresponds to the fraction of methyl methacrylate in the poly- (methyl methacry- 
late)(PMMA)/polystyrene(PS) diblocks used by Thurn-Albrect et al.^. This leaves three 
parameters whose effects can be explored. We have chosen to vary E^xt and do, and have 
examined two surface field strengths, AA^ = 0.2, and 0.5. To understand whether these 
surface fields are weak or strong, the surface energy of n polymers subject to these fields, 
Usu, can be compared to the electrostatic field energy of these polymers in volume Q. The 

10 



former of these energies is Usu = nkBTNA, while in the bulk cylindrical phase, the latter is 
Uei = KoeoE^^^n/2, where 

fiO = i^AfA + f^Bi^- Ia), (22) 

is the average dielectric constant of the system. Thus we consider the ratio 

Uel tQKQE^^p. 



Usu 2kBTANn 

-'ext 
2' 



2ANS 
where we have introduced the convenient scale of electric field 

and the value of the external field in these units, Eq = E^xt/^- At typical temperatures, 
T ^ 4307^, and for a typical volume per polymer chain of 100 nm^, this unit of electric 
field is £^ ~ 82V//im. At experimental temperatures around 160°C the dielectric constants 
appropriate to the PMMA-PS copolymer, with PMMA being the A block and PS the B 
block are^i^S^ ka = 6.0 and kb = 2.5, from which kq ^ 3.52. Thus Uei/Usu ^ 1.8Ei/AN. 
At external fields of order of tens of volts per micron, the surface field energies we consider 
are comparable to the electrostatic energies, and are therefore smaller than the rather large 
surface fielda^i of the experiments of Xu et al.— . 

We first consider the system with surface fields of AA^ = 0.5. We have calculated, for a 
given thickness, do, the excess free energies of various parallel, perpendicular, and interme- 
diate phases. That with the lowest excess free energy is the equilibrium one. Rather than 
plot the excess free energy, a, directly vs. applied field, we note from eq. El and the fact 
that the electrostatic contribution to the bulk free energy of the cylindrical phase is simply 
— eoKo-E'g^t/2, that for determining the globally stable phase, we can equally well plot the 
dimensionless surface free energy 

^"(^»)-;^ + r°^»- (^«' 

The derivative of this function with respect to the dimensionless external field Eq can be 
found from eq. d to be the negative of the dimensionless excess displacement field 



D.. 



dotoS 
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= -^^-™. (26) 

In fig 1. we show the excess free energy /„ for several phases as a function of the dimen- 
sionless external electric field for a thickness cLq = 7.5N^^'^a. First we note, as anticipated, 
that the surface free energy increases with external field, indicating that the surface contri- 
bution to the displacement is, in fact, negative i.e. in the opposite direction of the applied 
field. Secondly, we see that at this thickness, the stable morphology at small external fields 
is one with five layers of cylinders parallel to the plates. As the field is increased, a first order 
transition occurs to a mixed state in which the surface contribution to the displacement field 
is less negative. At sufficiently large fields, there is a second transition to the state in which 
all cylinders are parallel to the external field. The excess displacement field is even less 
negative, but is not zero due to the distortion of the cylinders near the plates. By repeat- 
ing such calculations as a function of thickness, we obtain the phase diagram shown in fig. 
2. Because the film thickness is a smooth function of the chemical potential, this diagram 
reflects much of what was anticipated earlier in our discussion of the phase diagram as a 
function of electric field and chemical potential. Note that for films of sufficient thickness to 
accommodate more than two layers of parallel cylinders, the surface field is sufficiently strong 
to impose a surface ordered layer which results in the existence of an intermediate phase. At 
higher fields, this surface ordered layer becomes too costly, and disappears resulting in the 
perpendicular phase. As anticipated, the electric field which brings about this latter transi- 
tion varies relatively weakly with film thickness. For sufficiently thin films, both orientations 
of cylinders, characteristic of the intermediate phase, cannot be maintained leading to the 
elimination of the intermediate phase. 

In determining the phase boundary between the perpendicular and intermediate phases, 
and between the perpendicular and the cylindrical phase with two parallel layers, we have 
set the characteristic length, e, of the surface field to a constant value 3N^^'^a/l6. It is not 
difficult to see from eqs. 1201 12] that if e is to equal an integer number n of grid points, of 
which there are 192 in the x direction, then the thickness do/N^'^a must take the discrete 
values 36/n. These discrete points are shown in fig. 2 where they are connected by a spline 
fit, shown dotted, to guide the eye. In contrast to the relatively smooth variation with 
thickness of the electric field at the transition from the perpendicular phase, that at the 
transition from the parallel phase varies rapidly with thickness. Hence the discrete values of 
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thickness imposed by a constant e is restrictive. Further, the behavior of this critical field 
with thickness has already been anticipated from our discussion of the thermodynamics of 
the transition. For these transitions, therefore, we chose e to be a constant fraction of film 
thickness, e/d^ = 1/32. It is readily seen, again from eqs. [201 1^ that this choice permits 
the film thickness do/N^^'^a to vary continuously. We have compared the two values of the 
critical fields obtained from these two different choices of e at the discrete values of do = 36/n 
permitted by a constant e, and found essentially no difference. 

We next present density profiles of the various phases for a distance between plates of 
do = l.haN^'"^. When the applied field is zero, the stable phase is that shown in fig. 3 
which we denote as having five parallel layers. The cut is in the x (vertical), y (horizontal) 
plane, a cut normal to the plane of the plates. The gray scale shows the volume fraction 
of A monomer from 0.0 (lightest) to 1.0 (darkest) in four bins of width 0.25. The A block 
is favored by the surface field of strength AA^ = 0.5. The distortion of the cylinders at the 
surface is clear. A cut in the yz plane, parallel to the plates, (not shown), confirms that 
these are, indeed cylinders with axis along the z direction. 

With an applied field of sufficient strength, the intermediate state becomes globally stable. 
A cut in the xy plane is shown in fig. 4 for a value Eo = 0.49. Most of the film thickness 
is occupied by cylinders which are perpendicular to the plates, and a cut in the yz plane 
halfway between the plates, (not shown), confirms that the cylinders are in a hexagonal array. 
It might be expected from fig. 4 that, adjacent to the surface, one would find a distorted 
layer of cylinders, just as in fig. 3. That this is not the case is shown by a cut through 
the intermediate phase in a yz plane very near the surface, fig. 5. One sees instead that 
there is an hexagonal array, such that the surface-favored A monomers occupy most of the 
surface layer. Clearly the hexagonal arrangement on the surface derives from the hexagonal, 
perpendicular, arrangement of cylinders in most of the film. Comparison of a sequence of 
yz cuts confirm that the symmetry axes of the upright cylinders in this intermediate phase 
coincide with the symmetry axes of the hexagonal array at the surface; i.e. the A-rich cores 
of the cylinders stand directly above the S-rich centers of the surface pattern. That the A 
monomers near the surface are located outside of the cylinders thereby occupying a larger 
area in the unit cell is clearly an effect due to the surface field. Because the A monomers are 
the minority component, Ja = 0.29, this configuration cannot extend too far into the film. 

For larger values of the electric field, the perpendicular phase is attained. Fig. 6 shows 
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the density profile in the xy plane at a field of Eq = 0.78. One clearly sees the distortion of 
the cyhnders produced by the surface. 

We next show, in fig. 7, the phase diagram for the weaker surface field of strength 
AA^ = 0.2. Again, we chose e = 3N^^'^a/16. One sees that the weaker surface field can 
no longer impose a surface morphology which differs from, yet coexists with, that of the 
remainder of the film, even for ones sufficiently thick to accommodate seven parallel layers. 
Thus there is no intermediate phase. Instead, with increasing external field, transitions from 
all cylindrical phases shown, up to seven parallel layers, are directly to the perpendicular 
phase. The maximum critical fields of the transition between parallel and perpendicular 
phases vary as l/v^O) as expected. This behavior is shown in the inset. The six maximum 
critical fields are shown by the open triangles which fall on a straight line in this plot. 

IV. DISCUSSION 

We have applied self-consistent field theory to a planar film of cylindrical-forming block 
copolymers. The plates between which the film is adsorbed prefer the minority A component. 
In the absence of an electric field, this causes the morphology to be that of cylinders whose 
axes are parallel to that of the plates. We have studied the equilibrium phase diagram in 
the presence of an applied electric field which tends to align the cylinders perpendicular to 
the plates, and have solved the Maxwell equations without approximation. The behavior as 
a function of field and film thickness is, for the most part, similar to that of lamellar-forming 
dibocks. As argued above, this is because the underlying physics, that of the competition 
between bulk and surface phases, is the same. In particular for weak surface fields or thin 
films, the surface field cannot impose a surface morphology which differs from the rest of 
the film, so that there is no intermediate phase. The first-order transition with increasing 
electric field from the phase of parallel cylinders to that of perpendicular cylinders is direct. 
However for stronger surface fields and thicker films, an intermediate phase does appear, one 
which is characterized by a boundary layer near the plates. 

The nature of this boundary layer is very different for cylinder-forming diblocks than for 
lamellar-forming ones for which this layer is simply a modulated lamellae. We considered 
different possibilities for this boundary layer. Because of the experimental report of a signal 
indicating the presence of both perpendicular and parallel cylinders^, we considered whether 
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a layer near the plates might consist of parallel cylinders, but found that it did not. In 
retrospect, this is not difficult to understand. First, the perpendicular cylinders aligned 
with the field could not connect in a simple periodic manner with cylinders aligned parallel 
to the plates as the two arrays would be incommensurate. Thus one of the two arrays would 
have to be distorted, with a concomitant increase of free energy. Second, the substrate does 
not favor parallel cylinders per se, rather it favors one monomer over the other, the minority 
A monomer in our case. In zero electric field, parallel cylinders are favored over perpendicular 
ones because the area of A monomers presented to the plates is larger in the former than in 

1 /2 

the latter. The ratio of areas is Aparf-Aperp ~ V/a > ^- In a non-zero field which favors 
the intermediate phase, however, the area presented to the plates by parallel cylinders is 
less than that presented in the configuration of hexagonal symmetry which we obtained and 

-1 /o 

which is shown in Fig. 5. In this case the ratio of areas is Aparf-Aint ~ /a /(I ~ Ia) < 1- 
We also note that this particular pattern of the surface layer derives from the fact that the 
minority A component is favored by the surface. Were the B component favored, as is the 
case in the system of Xu et aL^-^, one would expect the minority core of the cylinders to 
extend all the way to the plates over a range of surface fields, something which could be 
technologically useful. Only for very strong fields would the plates be essentially covered by 
the majority B component which would require the interior of the cylinders containing the 
A component to be truncated near the surface. This observation is in agreement with the 
dynamic self-consistent field theory results as seen in Fig. 1 1 (b) of Xu et al.— and Fig. 5 of 
Lyakhova et al.— . 

We have found a transition from the intermediate phase to the perpendicular phase for 
which the electric field needed to bring about the transition varies only weakly with film 
thickness. It would seem that this could be identified with the transition observed by Thurn- 
Albrecht et al^. There are, however, two difficulties with this interpretation. The first is that 
the electric field values at the transition are, in our calculation, larger than those observed 
in experiment. Further, the strength of the surface field in our calculation is weaker than 
that estimated in experiment. Were we to calculate critical electric fields for the stronger 
substrate fields of experiment, we would obtain even larger values as the critical electric 
field is expected to increase as the square root of the surface field. This comparison is 
reminiscent of that between calculate d^^i^^ and experimental^ electric field values for the 
transition from a body-centered cubic phase to a cylindrical phase in which the former were 
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also much larger than the latter. Explanations for this have been proposedi^^ which may 
resolve this difficulty. Even so, there remains the fact that, in the calculated intermediate 
phase, there are no cylinders in the parallel orientation. Therefore one would not observe 
in equilibrium the existence of both parallel and perpendicular cylinders along the phase 
coexistence of the intermediate and perpendicular phases. The resolution of this difficulty 
may be that, as the system cools from the disordered into the ordered phase, cylinders are 
nucleated near the plates in a configuration parallel to them while the major part of the ffim 
nucleates cylinders aligned with the electric field. Rearrangement of cylinders near the plates 
into that boundary layer which we predict to be the equilibrium one would presumably be 
a very slow process. If the boundary layer did not attain its equilibrium configuration, the 
experimental signal which detects both orientations of cylinders could be understood. 
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conversations. This work was supported by the U.S. -Israel Binational Science Foundation 
(B.S.F.) under grant No. 287/02 and the National Science Foundation under grant No. 
DMR-0503752. 



16 



1 T. L. Morkved, M. Lu, A. M. Urbas, E. E. Ehrichs, H. M. Jaeger, P. Mansky, and T. P. Russell, 
Science 273, 931 (1996). 

2 M. Park, C. Harrison, P. M. Chaikin, R. A. Register, and D. H. Adamson, Science 276, 1401 
(1997). 

3 S. Walheim, E. SchafFer, J. Mlynek, and U. Steiner, Science 283, 520 (1999). 

"^ T. Thurn-Albrecht, J. DeRouchey, and T. P. Russell, Macromolecules 33, 3250 (2000). 

^ T. Thurn-Albrecht, J. Schotter, G. A. Kastle, N. Emley, T. Shibauchi, L. Krusin-Elbaum, 

K. Guarini, C. T. Black, M. T. Tuominen, and T. P. Russeh, Science 290, 2126 (2000). 
^ K. Amundson, E. Helfand, X. Quan, S. D. Hudson, and S. D. Smith, Macromolecules 27, 6659 

(1994). 

T. Arnold, R. K. Thomas, M. Castro, S. M. Clarke, L. Messe, and A. Inaba, Phys. Chem. Chem. 

Phys. 4, 345 (2002). 
* K. Amundson, E. Helfand, X. Quan, and S. D. Smith, Macromolecules 26, 2698 (1993). 
^ A. V. Kyrylyuk, A. V. Zvelindovsky, G. J. A. Sevink, and J. G. E. M. Fraaije, Macromolecules 

35, 1473 (2002). 
^^ K. S. Lyakhova, A. V. Zvelindovsky, and G. J. A. Sevink, Slow Dynamics in Complex Systems, 

M. Tokuyama and I. Oppenheim eds.. Springer p. 217 (2004). 
" K. S. Lyakhova, A. V. Zvelindovsky, and G. J. A. Sevink, Macromolecules 39, 3024 (2006). 
12 G. G. Pereira and D. R. M. Williams, Macromolecules 32, 8115 (1999). 
1^ Y. Tsori and D. Andelman, Macromolecules 35, 5161 (2002). 
1^ R. Pandit, M. Schick, and M. Wortis, Phys. Rev. B 26, 5112 (1982). 
1^ M. W. Matsen, Phys. Rev. Lett. 95, 2583021 (2005). 
1^ T. Xu, A. V. Zvelindovsky, G. J. A. Sevink, K. S. Lyakhova, H. Jinnai, and T. P. Russell, 

Macromolecules 38, 10788 (2005). 
^■^ C.-Y. Lin, M. Schick, and D. Andelman, Macromolecules 38, 5766 (2005). 
IS M. W. Matsen, J. Chem. Phys. 124, 0749061 (2006). 

19 Y. Tsori, D. Andelman, C.-Y. Lin, and M. Schick, Macromolecules 39, 289 (2006). 
'^^ L. D. Landau, E. M. Lifshitz, and L. P. Pitaevskii, Electrodynamics of Continuous Media (Perg- 

amon Press, Oxford, 1984). 

17 



21 M. W. Matsen, J. Chem. Phys. 106, 7781 (1997). 

22 K. O. Rasmussen and G. Kalosakas, J. of Polymer Sciences: Part B: Polymer Physics 40, 1777 
(2002). 

23 G. Tzeremes, K. O. R. T. Lookman, and A. Saxena, Phys. Rev. E 65, 041806 (2002). 
2^ K. O. Rasmussen, J. of Polymer Sciences: Part B: Polymer Physics 42, 3695 (2004). 

25 T. Xu, A. V. Zvelindovsky, G. J. A. Sevink, O. Gang, B. Ocko, Y. Zhu, S. P. Gido, and T. P. 
Russell, Macromolecules 37, 6980 (2004). 

26 Y. Tsori, F. Tournilhac, D. Andelman, and L. Leibler, Phys. Rev. Lett. 90, 145504 (2003). 

27 P. Mansky, Y. Liu, E. Huang, T. P. Russeh, and C. Hawker, Science 275, 1458 (1997). 



18 



V. FIGURE CAPTIONS 

Fig. 1 Dimensionless excess surface free energy, /„, eq. 123 as a function of applied field, 
£^0) for thickness do = 7.5N^^'^a. The value of x^ = 18; and the surface field is AA^ = 0.5. 
The free energy of the perpendicular phase is shown with a solid line, that of the intermediate 
phase with a dashed-dotted line, and that of the parallel phase of five layers is shown with 
a dotted line. 

Fig. 2 Phase diagram as a function of electric field, Eq, and thickness, do/N^^'^a. The 
value of x^ = 18; ^"^^ the surface field is AA^ = 0.5. Parallel phases are denoted by a 
roman numeral corresponding to the number of cylinders in the film. The intermediate and 
perpendicular phases are marked "Intermed", and "Perp" respectively. Discrete calculated 
points on the boundary between perpendicular and intermediate phases are shown by solid 
dots. The dotted line is a spline fit to them. 

Fig. 3 Density profile in the xy plane of the parallel phase with five layers of cylinders. 
The thickness is do/N^'^a = 7.5, and the applied field is zero. The gray scale shows the 
volume fraction of A monomer from 0.0 (lightest)to 1.0 (darkest) in four bins of width 0.25. 

Fig. 4 Density profile in the xy plane of the intermediate phase. The thickness is 
do/N^^'^a = 7.5 and the applied field is now Eq = 0.49 

Fig. 5 Density profile of the intermediate phase shown here in a yz plane very close to 
the plates. 

Fig. 6 Density profile of the perpendicular phase in the xy plane. The thickness is 
do/N^^'^a = 7.5 and the applied field is now Eq = 0.78 

Fig. 7 Phase diagram as a function of electric field, Eq and thickness d^/N'^/'^a for a 
surface field AA^ = 0.2 and x^ =18. Inset shows as open triangles six calculated maximum 
electric field values at the transition from parallel to perpendicular phases plotted vs. I/c/q . 
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Chin- Yet Lin and M. Schick, Fig. 3 
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Chin- Yet Lin and M. Schick, Fig. 4 
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Chin- Yet Lin and M. Schick, Fig. 5 
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Chin- Yet Lin and M. Schick, Fig. 6 
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Chin- Yet Lin and M. Schick, Fig. 7 
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